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Abstract 



U 

<. 

, The zero temperature in plane penetration depth in an untwinned single crystal of optimally 

doped YBa2Cu306.93 is highly anisotropic. This fact has been interpreted as evidence that 
a large amount of the condensate resides in the chains. On the other hand, the temperature 
dependence of Aq(T)/Aq(0) and Ab(T)/Ah(0) (where the ^-direction is along the chains) are not 
very different. This constrains theories and is, in particular, difficult to understand within a 
' proximity model with d-wave pairing only in the Cu02 plane and none on the CuO chains but 

instead supports a more three dimensional models with interplane interactions. 

o , 

^ ■ Introduction 

' Most theories of the high Tc oxides start with a Cu02 plane which is the common building block 

^ ■ found in all the copper oxide superconductors. Such a system is tetragonal and is sometimes 

', described by simple two dimensional tight binding bands with first and second nearest neighbour 

' hopping. The mechanism involved in the superconductivity is still unknown but there is now strong 

O ■ evidence, if not yet a consensus, that the gap has d-symmetry|, |, |, |, |, |, Q Of course there 

are other elements to the structure of the typical copper oxide superconductor but a stack of Cu02 
planes weakly coupled through a transverse hopping can be taken as a simplified first model. 
^ • While in many of the oxides t± is small — perhaps of the order of 0.1 meV in Bi2Sr2Cu Cu208,|ql 

and of the order of a few meV in LaSrCu04 — in YBa2Cu307_5 (YBCO) at optimum doping, it is 
much larger and of the order of a few tens of meV which is almost of the same order of magnitude 
as the in-plane hopping integrals and indicates that this material may be fairly three dimensional. 
In addition, YBCO has chains (CuO) as well as planes (Cu02). In such circumstances, the system 
can be expected to be significantly orthorhombic with the source of orthorhombicity residing in the 
chains. In an orthorhombic system, the irreducible representation of the point group crystal lattice 
which contains the dx2_y2 part also contains s^2_^_y2 and So (constant) parts, and these can mix in 
the gap so that we cannot expect a pure d-wave order parameter. 

Strong evidence that the chains play a very important role in optimally doped YBCO is obtained 
from infrared and microwave experiments on untwinned single crystals of YBa2Cu306.93. Far 
infrared experiments!^ can be used to measure the absolute value of the in plane penetration 
depth at zero temperature in each of the two principle directions denoted by a and b with the 
chains oriented along b. The results are Afe ~ 1030A and Xa ~ I6OOA for YBa2Cu306.93 with 
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Figure 1: Energy dispersions along the line ky — for = 0. There is an avoided crossing at = ^2- 
The band parameters are typical and are ti = 100 nicV, <2 = 80 meV, /ii = —80 meV, /i2 — 40 meV and 
t_L = 25 meV. 



Tc =93. 2K. On the other hand, the temperature dependence |1C, 0] of the normahzed penetration 
depth, A(r)/A(0), is almost the same in both directions. These results have been taken as evidence 
that the chains carry a significant amount of the condensate density and that the gap in the 
chains is of the same order of magnitude as in the planes, a fact confirmed in current-imaging 
tunneling spectroscopy (cits) experiments. [14, 15] It would also follow from the observed large 
orthorhombicity that the order parameter will not be of a pure (1^2 _y2 symmetry as previously 
stated and repeated here. It should contain a significant s-admixture due to the existence of the 
chains which are coupled to the planes and participate importantly in the superconductivity. 



Proximity Model 

A first model, which can be used to get some insight into the situation for optimally doped YBCO, is 
that of planes and chains coupled through a transverse tunnelling matrix element t± with the pairing 
interaction assumed to reside exclusively in the CO2 plane. The simplest electronic dispersion 
relations for such a model |jl^, ^ |l^ are: 

^1 = -2ti[cos{kx) + cos{ky)] - ^ii (1) 
6 = -2t2Cos{ky) - fj.2, (2) 

where ti is the first neighbour hopping in the Cu02 plane and t2 in the CuO chains with ^ the 
chemical potential. Application of an interplane matrix element t± will mix the plane and chain 
band with resulting band energies having the form: 



.± = ^±^(^)%t,(k)^ (3) 

with 

t±{k) =2t±cos{k^/2) (4) 
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Figure 2: The Fermi surface is shown for a range of kz between fcz = and kz = tt. When kz = tt, the 
chain-plane couphng vanishes and the two pieces of Fermi surface are those of the isolated chains and planes. 
As the chain plane coupling increases, the Fermi surfaces hybridize and are pushed apart. The amount of 
hybridization at a given k depends on the relative sizes of and (^i — ^2)^- There is an avoided crossing of 
the Fermi surfaces when ^1 = ^2 = 0. 



in these relationships kx, ky and kz are components of momentum ranging from (— 7r,7r) in units 
of the inverse of the lattice parameters. Results are shown in Fig. 1. As seen, the effect of the 
plane-chain coupling is largest where ^i(k) = ^2(k) at which point there is an avoided crossing. 
In Fig. 2, wc show how the Fermi surfaces arc pushed apart in A;-space and sec that the amount 
of distortion of the Fermi surfaces depends on their proximity to the avoided crossing in the 2-D 
Brillouin zone. The various contours are for different values of kz with kz = tt corresponding to no 
coupling {t±{kz) = 0) and A;^ = corresponding to the largest effect. The area between these two 
outer contours corresponds to the Fermi surface dispersion in the z-direction. If it was a perfectly 
flat cylinder-like structure in this direction, it would project into a single curve in the 2-D Brillouin 
zone. 

One result of the large critical temperature in the oxides is that they have an extremely small 
coherence length ^o- The coherence length is the distance scale over which the superconducting 
order parameter A may vary spatially. In the BCS theory it is : 

where Vf is the electron Fermi velocity, and A is the magnitude of the BCS order parameter. For a 
BCS superconductor with = 100 K, A = 15 meV. The Fermi velocity can be estimated from the 
bandwidth of the conduction band, and since the high Tc materials are highly anisotropic, there will 
be substantial differences between Vf m. the various directions. The coherence length will therefore 
be anisotropic as well. In the a and h directions (parallel to the Cu02 planes) the Cu02 bandwidth 
is t\\ ~ 1 eV. The Fermi velocity can be estimated as: 

at\\ 

Vfa ^Vfb^ ^ 
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Figure 3: The quasiparticle excitation spectra E± are plotted along the line kx = (solid curves). For these 
curves, t± = 20 meV. The excitation spectra in the decoupled {t± = 0) limit are plotted for comparison 
(dashed curves). 



where a is the lattice constant in the a-direction. For the BSCCO compounds, a 
coherence length in the planes is therefore: 



SA. The 



atii 



TT 



2 A 



25A. 



The bandwidth along the c axis is considerably smaller than in the Cu02 planes. In the BSCCO 
compounds, it is typically t± ~ 10 meV. Since the unit cell size along the c-axis is c ~ 30A, the 
coherence length in BSCCO is 



ct± 



2A. 



These values are typical for most of the high Tc materials and distinguish them from the conventional 
materials in one important way: the fact that is substantially less than the length of the unit cell 
along the c axis allows the order parameter to vary spatially over the unit cell. In the conventional 
materials, where the coherence lengths are 10^ or lO^A, the structure of the unit cell is invisible 
to the gap. In the high Tc materials, the value of the gap may depend on the layer type: chain or 
plane. 

Next we wish to include the superconductivity. We will assume that the pairing is operative 
only in the Cu02 plane and that the chains become superconducting only through the tunnelling 
matrix element t±. In mean field theory, the Hamiltonian, H, is: 



H-Nii = Y, C+(k)Q(k)C(k) + const.. 



(5) 



with: 



C{k) 



cikT 

C2kT 



(6) 
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Figure 4: The critical temperature in the intrinsically superconducting planes is lowered by the proximity 
of the intrinsically normal chains. 



where c\■^^^ are the creation operators for electrons in layer i and spin a. Here the matrix Q(k) is 
4x4 and equal to: 

"Ci(k) -Ak m 



Q(k) 



-6(-k) 







The gap Ak is given by: 



t*(k) 6(k) 
-t(-k) -6(-k) 

= ^ X! H;k'(ci-k'iClk'T)! 



(7) 



(8) 



where Q is the volume, Vkk' is the pairing potential in the copper oxide plane denoted by the 
subscript 1. We will assume Vkk' to be separable and have d-wave symmetry. That is: 



with: 



77k = cos(A;^) - cos(A;2;). 



(9) 



(10) 



Diagonalisation of the Hamiltonian (5) leads to four energy bands Ei = E^,E2 = E-,Es = 
-E^,E4 = -E+ with: 



(i'-a6)'-(6Ak)2. 



(11) 



In Fig. 3, wc show the energy bands for t± = (dashed curves) and t± = 20 meV (solid curves). 
For the uncoupled situation shown for comparison (dashed curve), the chain band exhibits no gap 
while the plane band shows the usual BCS gap at its Fermi surface. Note that the curves are for 
fcj; = as a function of ky and that the d^2_y2 gap of the form (10) is nonzero along that line in the 
two dimensional Brillouin zone. When t± ^ 0, the bands mix as they would in the normal state 
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Figure 5: The value of the lower energy band E^{k) (dashed line) is plotted along the Fermi surface (solid 
line) of the chain-plane model for = 0. The chain-plane coupling strength is t±_ — 25 meV. 



with a band crossing. In addition, a gap is induced in the E- branch corresponding to induced 
superconductivity in the chains coming from the proximity matrix element t^. This gap is small 
compared with the value of the original gap in the Cu02 plane which is seen at somewhat higher 
momentum in this figure. 

A first question that needs to be answered is how is Tc changed when t±_ is switched on. This 
is shown in Fig. 4 where we have plotted the value of as a function of t± for a case when the 
critical temperature Tc =100 K in the limit t±_ = 0. We see that as t±_ (in meV) increases, Tc is 
reduced by the proximity of the chains but remains substantial even for a value of t±_ = 20 meV. 

In Fig. 5, we show the chain-plane Fermi surfaces (solid curves) as contours in the first Brillouin 
zone for fc^ = and t±_ = 25 meV. The value of E-{k) (dashed line) is also plotted along the Fermi 
curves. It is seen that in the region of the chain Fermi surface which is not close to the plane Fermi 
surface, i.e. the region around (vr, 7r/2) in the figure, the induced gap is small. In fact, it can be 
shown that in this region: 

[ei(k)-e2(k)]2»t^(k)2, (12) 

and that: 

^_(k) ~ r ' 



' '"(6-6)2 

~ 0.1 meV, (13) 

for the parameters used in the calculations. This small gap is an important generic feature of a 
pure proximity model with a plane-chain Fermi surface (2 sheets). 

In a coupled two band model, the expression for the electromagnetic response tensor K^^ is 



complicated and has the form [17|: 



where: 

e2 



K^.u = ^ {gV(0,0;0)|^^^ - G^R^Mm) , (14) 



^^(O'OSO) = -7T7^ E E[%(k,k)].,[7^(k,k)b, 
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Figure 6: In-plane penetration depth for a d-wave order parameter, the penetration depth in the a direction 
(perpendicular to the chains) is nearly that of a pure d wave superconductor in the absence of chains. The 
penetration depth in the b direction has a very different temperature dependence from that in the a direction 
because the size of the induced gap in the chains is much different from the size of the gap in the planes. The 
relative bandwidths of the chains and planes were determined by setting A^(0)A^(0) ^ 2.5, in accordance 
with experiment. 



dEi 



J-{Ei)-f{Ej 



E, 



(15) 



where e is the charge on the electrons, c is the speed of light, is the volume, h Planck's constant, 
and [7^(k, k)]jj is the appropriate electromagnetic vertex which, in the familiar one band model, 
would be the Fermi velocity. Here we have [7^(k, k + q)]jj = [/'l"(k)[7^(k, k + q)]jjC/(k + q) with 
U(k) the 4x4 unitary matrix that diagonalizes the Hamiltonian (|^. The vertex [7^(k, k)]jj is 
related to the dispersion curves in the bands: 



[7M(k,k)] 



-l)^-iVkQ°(k) 



(16) 



where Q*^ is the Hamiltonian matrix (Eq. (|^)) in the normal state. In the limit of a single band 
( [TsD properly reduces to the familiar result: 



K. 



e 



dEk 9ek 



(17) 



where v(k) is the Fermi velocity and E]^ = [e^ + A^]^/^. The second term in equation (^) is 
evaluated in the normal state and is related to the value of the zero temperature penetration 
depth. 

Results for and A^'^ coming from the chains and planes in our model are shown in Fig. 6. In 
our notation, the chains are along the b-direction. For the current along a-direction, the penetration 
depth (solid curves) shows linear dependence at low temperature as expected for d-wave. This is to 
be contrasted with the exponentially activated behaviour found for the s-wave case. For the currents 
in the b-direction (dashed curve), the chains contribute significantly to the superfluid density and 
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Figure 7: A comparison of penetration depths found by experiment and by the proximity model. The 
experimental data |12| on untwinned crystals of YBCO is shown in (a). Squares are Aa(0)^/Aa(T)^ and 
triangles are AQ(0)^Af,(r)^. In (b), the theoretical values of the penetration depth are shown for the 
proximity model in the a (solid curve) and b (dashed curve) directions. 



the shape of the curve is very different. It shows a strong upward bending at low T which reflects 
the low energy scale noted in formula ([l^ for the value of the gap which comes from regions of the 
chain Fermi surface well away from the band crossing point. This feature is generic to proximity 



models with pairing confined to the Cu02 plane. The data of Bonn and Hardy pO|, 11, 12] in pure 
single untwinned crystals of YBCO are shown in Fig. 7 where they are compared with our results. It 
is clear that no low energy scale is observed in the data along the b-direction. In fact, on normalized 
plots, the observed temperature variation is very similar between a- and b-directions with the slope 
in the b-direction slightly steeper. This indicates clearly that the gap in the chain is large and that 
no small energy scale exists. To understand better how the low temperature slope is related to the 
d-wave gap, we show, in Fig. 8, results for the penetration depth in a one band model with gap 
^Ao[cos(A;a;) — cos{ky)] for several values of the ratio 2Ao/A;_b Tc. In obtaining these results, a BCS 
temperature variation was taken for the temperature dependence of the gap. It is clear that small 
gap values give a curve with concave upward curvature while for large values the curve is concave 
downward. For BCS 2Aq/ ksTc = 4.4 and the curve is nearly a straight line with 45° slope. If for 
the same zero temperature gap value the critical temperature is decreased, the curve will clearly 
have a smaller slope. This is because in that case the system at low temperature is expecting that 
it should have a straight line behaviour with an intercept at much higher temperature. 

An analytic result that can be proved for a simple d-wave model with circular Fermi surface in 
the 2-dimensional Brillouin zone and gap variation of the form 



A = Ao COs(2(/)) 



(18) 
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Figure 8: The penetration depth, [A(0)/A(r)]^, for different values of 2zi,„ax/7c. From top to bottom the 
values are: 8, 7, 6, 5, 4.4, 3 and 2. The data points are from microwave experiments. [l^ 



where (p is an angle along the Fermi surface is that 

47rne2 21n(2) 



1 



(19) 



for T — 5- 0. This shows that the slope as a function of reduced temperature is inversely proportional 
to the ratio 2Ao/kBT. As this ratio increases, the slope becomes less steep. Here n is electron 
density and m is electron mass. 



Interband Pairing 

If one wants to remain within a model where there is no pairing interaction in the chains, one way 
to increase the value of the chain gap is to include off diagonal pairing in a two band BCS model. 
The BCS equations in this case are [18-19] 

1 



= X! (^k,q,llXq,l + Vk,q,12Xq,2) 



^k,2 = X! (^k,q,12Xq,l + ^k,q,22Xq,2), 



(20) 



where the pairing susceptibility is: 

Xq,c 



(CqT. 



o'--qJ.,Q/ 



2^q,a ^ \2kBTj 



q,a 



tanh 



(21) 



for the a'th band with the brackets (•) indicating a thermal average of the pair annihilation opera- 
tors. In calculations, we will assume V22 to be zero but take V12 = V21 finite. This is the term that 
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Figure 9: Model of YBCO Fermi surfaces for chains (long dashed curve) and planes (closed short dashed 
curve) in the first Brillouin zone. The (tt, tt) point is at the center of the figure. For the chains the parameters 
{ta, Ca, Ba, fj,a} in Eq. H are {-50, -0.9, 0, 1.2} and for the planes they are {100, 0, 0.45, 0.51}. 



couples the chains and planes and makes the chains have a gap. For the electron dispersions, we 
take a model vi^ith up to second neighbour hopping with: 



Ck,a = -2ia [(1 + ea) cos{kx) + cos(A;y) 

— 2Ba cos{kx) cos{ky) — (2 — 2Ba — fJ-a)] 



(22) 



where the two new parameters not in equation (1) and (2) are the second neighbour hopping Bq, 
and the orthorhombic distortion e^. As a model we take {ta, ^a, Ba, fJ-a} to be {100,0,0.45,0.51} 
for the planes and {—50, —0.9,0, 1.2} for the chains. The resulting Fermi surfaces for chains (long 
dashes) and planes (short dashes) are shown in Fig. 9. 

To solve for the gaps of equations (^0[) , we need to make some choices for the pairing potential 
Vk,k',a/3- In the previous section, we chose a separable form. Here we use a different alternative, 
which is based on the nearly antiferromagnetic Fermi liquid approach, and assume: 



Vk,q,a/3 = gafS 



-t 



l+^2|k_q_Q| 



(23) 



which is proportional to the antiferromagnetic spin susceptibility with magnetic coherence length ^ 
taken from Millis, Monien and Pines [pi] ] and t = lOOmeV sets the energy scale. In equation (p3|), 
Q is the commensurate wave vector (vr, vr) and therefore the repulsive interaction (^) is peaked 
at (vr, vr). This interaction leads directly to a dx2_y2 gap in a single plane. If the coupling ga/s for 
a ^ P is different from zero superconductivity is induced in the chains and the gap no longer has 
pure dr^2_y2 symmetry in both chains and planes. It will be an admixture of: 



{d.^2 2) 



1 

cos{kx) + cos{ky) 
cos{kx) — cos{ky). 
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Figure 10: The order parameter in the first BriUiouin zone for (a) the plane layers and (b) the chain layers. 
Beside are the projections of the d components, (c) and (d), and the s components, (e) and (f). The vertical 
scale in all frames is the same. Note that the relative phase of the d-components in the two layers are the 
same while that of the s-components are opposite; this is caused by the interlayer interaction, 14c, q, 12, being 
negative (ie, repulsive). 



In Fig. 10, we show results for a case {511, 5121 922} = {26.2, 10, 0} with the critical temperature 
taken to be 100 K representative of the copper oxides. What is shown in frames (a) and (b) are the 
gaps as a function of momentum in the first Brillouin zone for the plane and chain, respectively. 
The projection on d-wave and s-wave manifold are also shown in (c) and (e) for the plane and (d) 
and (f) for the chain. We see that the orthorhombic chains can lead to a significant mixture of s 
and d components even for the plane case. 

A useful representation of these gap results is to show the contours of gap zeros on the same 
plot as the Fermi surface. This is presented in the series of frames shown in Fig. |ll[ The top 
frames apply to the planes while the bottom frames apply to the chains. In all cases, (a), (c), 
(e) for the planes and (b), (d), (f) for the chains, the same Fermi surface (dashed curves) was 
used. By choice, the Fermi contour have tetragonal symmetry in the top figure while the chain 
Fermi surface is a quasi straight line along as is expected for chains along y in configuration 
space. The pictures are for three different values of pairing potential. The first set of two frames 
(a) and (b) are for {511,^12,522} = {29.9,5,0}, i.e. very little coupling between chains and planes 
(off diagonal 512 small). In this case, the gap in the plane is nearly pure d-wave as is also the 
induced gap in the chains. As the coupling 512 is increased to {511,512,522} = {26.2,10,0} a 
significant s-wave component gets mixed into both solutions and the gap nodes move off the main 
diagonals of the Brillouin zone. (This is the solution that is plotted in Fig. 10). The gap nodes 
still cross the Fermi surfaces in both chains and planes. As the coupling is further increased to 
{51I) 5i2i 522} = {9.18,20,0}, Fig. 11 (e) and (d), the gap nodes move far off the diagonal and for 
the chains they no longer cross the Fermi surface so that there is a finite minimum value of the gap 
on this sheet of the Fermi surface. 

The amount of admixture of each of the three components in equation (^) are shown as a 
function of off diagonal 512 in Fig. 12 for the last two cases, namely 522 = 26.2 left frames and 
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Figure 11: The Fermi surface (dashed curves) and gap nodes (soUd curves) for a Cu02 plane layer (top 
frames) and a CuO chain layer (bottom frames) for three different interlayer interaction strengths (left, 
middle and right frames). As the strength of the interlayer interaction, gi2, is increased (left to right frames) 
the proportion of the s-component in both layers increases. If the interlayer interaction were further increased 
the gap nodes would leave the Brillouin zone altogether and the order parameter would become s-like. Note 
that the Fermi surface in the Cu02 plane layer is tetragonal but that the gap node is orthorhombic. 



gii = 9.18 right frames. The amplitude of the gap component involved is given in meV. Frames (a) 
and (c) are for the chains while frames (b) and (d) are for the planes. In all cases, the solid curve 
is the (i2,2_y2-component, short dashed the Sa^.2_|_y2-component and long dashed the So-component. 

For the first choice of intralayer interaction (left frames), gu = 26.2, and there is no order 
parameter in the chains when there is no interlayer interaction (ie, gi2 = 0) and the order parameter 
in the planes is pure d-wave. As the interlayer interaction is increased from zero, s-wave components 
appear in the planes and all three components appear in the chains. This "s + d mixing" is caused 
by the breaking of the tetragonal symmetry upon the introduction of the chains; there is no relative 
phase between the s- and d-wave components within either the planes or chains but there can be a 
relative phase between the order parameter in the planes and chains. In the range of gi2 explored 
here the d-wave component in the plane remains dominant but for sufficiently strong interlayer 
interaction the isotropic s-wave component eventually dominates |]l3| (ie, the gap nodes disappear). 
For interaction parameters {511,512,522} = {26.2,10,0} the critical temperature is lOOK and the 
maximum value of the gap in the Brillouin zone is 27.5meV in the planes and S.OmeV in the chains, 
while the maximum values on the Fermi surfaces are approximately 22meV and 7meV respectively. 
The ratio 2A max/Tc is 6.4 in the planes and 1.9 in the chains. 

For the second choice of intralayer interaction (right frames), gu = 9.18, there is no order 
parameter in either the chains or the planes when there is no interlayer interaction (ie, 512 = 0). As 
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Figure 12: Calculation of the zero temperature order parameters as a function of the interlayer interaction, 
gi2, for two fixed values of the interlayer interaction, gn, (left and right pairs of frames) presented for the 
planes, (b) and (d), and chains, (a) and (c), separately. In all frames the solid curve is the d-wave component 
of the order parameter, the short dashed curve is the extended s-wave component and the long dashed curve 
is the isotropic s-wave component. In the left frames gu = 26.2 and for gi2 — 10, Tc — lOQK. At gi2 = 
the order parameter is zero in the chains and is pure d-wave in the planes. As the interlayer interaction is 
increased the order parameter becomes present in the chains and there is a mixing of s-wave components. 
In the right frames gu — 9.18 and for 1712 — 20, — lOOK. At gi2 — the order parameter is zero in both 
the chains and the planes. As the interlayer interaction is increased the order parameter becomes present in 
both the chains and planes and there is a mixing of s-wave components with the isotropic s-wave component 
eventually becoming dominant. The feature at gi2 ~ 25 occurs when the gap node leaves the Brillouin zone. 
As discussed in the text there is a gi2 —gi2 symmetry. Both gu and gi2 are in units of ti; g22, the 
coupling in the chains, is set equal to zero. 



the interlayer interaction is increased d-wave and then s-wave components of the order parameter 
appear in both the planes and chains. Again, there is no relative phase between the s- and d-wave 
components within either the planes or chains but there can be an overall relative phase between 
the order parameter in the planes and chains. At approximately gi2 = 15 the gap nodes no 
longer cross the Fermi surface in the chains; the feature at gi2 ~ 25 coincides with the gap nodes 
leaving the Brillouin zone and isotropic s-wave becoming dominant. For interaction parameters 
{511)912)522} = {9.18,20,0} the critical temperature is again lOOK and the maximum value of 
the gap in the Brillouin zone is now 32.8meV in the planes and 20.1meV in the chains, while the 
maximum values on the Fermi surfaces are approximately 27meV and 17meV respectively. The 
ratio 2Z\niax/7c is 7.6 in the planes and 4.7 in the chains. 

Note that for gi2 > all of the s-wave components of the order parameters in both the planes 
and chains have the same relative sign and the d-wave components have opposite signs, while 
for gi2 < all of the relative signs are reversed but that the magnitudes of the components are 
insensitive to the sign of gi2- 

In Fig. |l^ the magnitude of the gap is plotted as a function of angle, 9, along the Fermi surface 
measured from the vertical. In frame (a) of Fig. ^ the local maxima of the gap on the Fermi 
surface are 16 and ISmeV; in (b) they are 1 and 7meV, and in (c) they are 25 and 3meV. In (d) 
one can see that there are no gap nodes which cross the Fermi surface; the maximum and minimum 
value of the gap on the Fermi surface are 17 and 4meV respectively. 

In Fig. |l^ we have plotted the magnetic penetration depth calculated with the lowest three 
harmonics (|2^) of the solutions of the BCS equations ( pO[ ) for the two choices of interaction pa- 
rameters. The solid curve is for the y-direction (along the chains) and the dashed curve is for 
the x-direction (perpendicular to the chains). The dotted curve is 1 — (T/Tc)^ and is plotted for 
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Figure 13: The magnitude of the gap on the Fermi surface as a function of angle for the four cases of 
Figs. The angle 9 is measured from the center or (tt, tt) point of the Brillouin zone with the y-axis (ie, 
the vertical) corresponding to = 0. Frames (b) and (d) do not span all angles due to the Fermi surface not 
being closed in the chain layer. For the first choice of interaction parameters, {ffii, 312, 522} ~ {26.2, 10,0}, 
the ratio 2Z\,„ax(FS)/21::j where Z\„-iax(FS) is the maximum value of the gap on the Fermi surface, is 4.3 and 
1.6 for the planes and chain respectively; for the second, {311,1712,522} = {9.18,20,0}, they are 5.7 and 3.8. 



comparison. The ratio \yy/\xx at zero temperature is 1.37 for both interaction parameter choices 
since the zero temperature penetration depth is a normal state property. The zero temperature 
penetration depth is largely governed by the bandwidth (ie, 4tQ,(2 — Ca)) - the larger the bandwidth 
the larger the zero temperature penetration depth. 

As pointed out above, the curvature of the penetration depth curve, X[^{T), is largely governed 
by the ratio 2Z\max/7'c and is a straight line for the d-wave BCS value of 4.4. The presence of the 
chain layer and the interlayer interaction increases this ratio in the plane layer but it remains low 
in the chain layer due to the absence of an interaction in this layer. It is this lower value that makes 
Xyy{T) (along the chains) have upward curvature (solid curves). Including pairing in the chains 
will push the solid curve towards the dashed one and make the initial low temperature slopes fall 
closer to each other. 



Conclusions 

In conclusion, the present data on single crystal untwinned YBCO at optimum doping suggest 
that the proximity effect incorporated into a single perpendicular tunnelling parameter t± cannot 
account for the observation. If interplane pairing is included the situation is greatly improved 
provided the off diagonal pairing is increased sufficiently to produce a gap on the chain which is of 
the same order of magnitude as that in the planes. Similar values of the gaps on the chains and 
planes can be taken as evidence that optimally doped YBCO is fairly three dimensional and that 
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Figure 14: The magnetic penetration depth for the two sets of interaction parameters. Which show the 
magnetic penetration depth along (solid curve) and perpendicular (dashed curve) to the CuO chains. The 
dotted curve is 1 — (T/Tc)^ and is shown for comparison. The chains, due to their Fermi surface, do not 
contribute appreciably to the penetration depth perpendicular to the chains (dashed curves). The ratio 
i^yy/^xx)"^ is 1.37 for both sets of parameters since this is a normal state property. 



the coherence length in the z-direction may not be sufficiently short to allow spatial inhomogeneities 
to exist within a unit cell. 
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